A=[1,2,-1;1,4,1;1,1,0;4,0,1];
b=[2;4;3;6];

X=zeros(3,8);
c=[3,-2,5];

z=zeros(1,8);

X(:,1)=[0;0;0];

X(:,2)=[0;0;1];

X(:,3)=[0;1;0];

X(:,4)=[0;1;1];

X(:,5)=[1;0;0];

X(:,6)=[1;0;1];

X(:,7)=[1;1;0];

X(:,8)=[1;1;1];

z=c*X; 

while (1)

    [m n]=find(z==max(z));

    if A*X(:,n)>=b

        z(n)=0;

    else

        x=X(:,n)

        zmax=z(n)

        break;

    end

end